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We present a statistical-model description of fission, in the framework of compound-nucleus decay, which 
is found to simultaneously reproduce data from both heavy-ion-induced fusion reactions and proton-induced 
spallation reactions at around 1 GeV. For the spallation reactions, the initial compound-nucleus population is 
predicted by the Liege Intranuclear Cascade Model. We are able to reproduce experimental fission probabili- 
ties and fission-fragment mass distributions in both reactions types with the same parameter sets. However, no 
unique parameter set was obtained for the fission probability. The introduction of fission transients can be off- 
set by an increase of the ratio of level-density parameters for the saddle-point and ground-state configurations. 
Changes to the finite-range fission barriers could be offset by a scaling of the Bohr- Wheeler decay width as pre- 
dicted by Kramers. The parameter sets presented allow accurate prediction of fission probabilities for excitation 
energies up to 300 MeV and spins up to 60 h. 
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I. INTRODUCTION 

Although seventy years have passed since the seminal 
works of Bohr and Wheeler (ill and Weisskopf and Ewing 
121 and the establishment of a qualitative understanding of the 
de-excitation mechanism of excited nuclei, quantitatively ac- 
curate and universally applicable models do not yet exist. This 
is partly due to the vastness of the amount of nuclear data that 
must be fed into the models, and partly to the uncertainties 
in the fundamental ingredients of de-excitation, such as level 
densities and emission barriers. Even the choice of the math- 
ematical formalism, however, is not devoid of confusion, as it 
was already pointed out by Moretto yll and Swiatecki |J4tl. 

One way to lift the degeneracy of the ingredients of the 
model is to explore diverse regions of the compound-nucleus 
parameter space. A systematic study of nuclei with different 
masses, excitation energies, spins and isospins would be sen- 
sitive to most of the assumptions of the de-excitation model. 
The long-term goal of such an investigation would be to iden- 
tify a minimal set of physical ingredients necessary for a uni- 
fied quantitative description of nuclear de-excitation chains. 

The production of excited compound nuclei can proceed 
from several entrance reactions. There has been a long history 
of compound-nucleus studies using heavy-ion-induced fusion 
reactions. These reactions allow one to specify the compound- 
nucleus mass, charge and excitation energy; however, a dis- 
tribution of compound-nucleus spins is obtained. Statistical- 
model parameters such as fission barriers are quite sensitive 
to spin. Heavy-ion-induced fission probabilities, evapora- 
tion spectra, residue masses can generally be reproduced in 
statistical-model calculations. However, some fine tuning of 
the statistical-model parameters to the mass region or reaction 
is often needed. A unified description over all mass regions 
is still lacking even for these reactions. Some work towards 



this, concentrating of the parameters describing the shape of 
the evaporation spectra, is presented in Ref. JSJl. 

Another typical entrance channel for the production of ex- 
cited compound nuclei is spallation. The present interest in 
spallation derives mainly from the applications to Accelerator- 
Driven Systems (ADS), namely accelerator-based reactors for 
the transmutation of nuclear waste. At incident energies rel- 
evant for transmutation (a few hundred MeV), an appropri- 
ate theoretical tool for the description of proton-nucleus reac- 
tions is the coupling of an intranuclear-cascade model (INC) 
with a nuclear de-excitation model. It is assumed in the INC 
framework that the incoming particle starts an avalanche of 
binary collisions with and between the target nucleons. When 
the cascade stage ends, an excited and thermalized remnant is 
formed, with a basically unchanged density. In the subsequent 
de-excitation stage, the remnant gets rid of the excess energy 
by particle evaporation and/or fission. For these reactions, 
the need for a model to predict the initial compound-nucleus 
mass, charge, excitation, and spin distributions adds some un- 
certainty in our ability to constrain the statistical-model pa- 
rameters by fitting data. However, spallation reactions allow 
us to explore different regions of compound-nucleus spin and 
excitation energy than can be probed with fusion reactions 
alone and thus can be important in parameter fitting. 

The role of a transient fission width is currently of some 
controversy. Fission transients are where the fission decay 
width is not constant, but increases from zero towards its equi- 
librium value 1I2I. Fission transients were first introduced to 
help explain the large number of neutrons emitted from a fis- 
sioning system before the scission point was reached 131 ■ The 
statistical model assumes there is an equilibrium is all degrees 
of freedom including the deformation degrees of freedom as- 
sociated with fission. If all compound nuclei have spherical 
shapes initially, then they cannot instantaneously fission as it 
takes as finite time to diffuse towards the saddle and subse- 



quently the scission point. The transient time, the time scale 
necessary for the system to explore large fluctuations in the 
deformation degrees of freedom, is a function of the viscosity 
of the nucleus. The predicted fission probability is also very 
sensitive to the assumed initial deformation 180 which may 
depend on the entrance channel. 

The transient time is often called a fission delay as fission is 
suppressed during this period. If the excitation energy of the 
compound nucleus is large enough, then there will be a prob- 
ability of emitting a light or possibility even an intermediate- 
mass fragment during the fission delay. Neutron emission 
lowers the excitation energy and charged-particle emission 
also lowers the fissility of the nucleus by increasing its fis- 
sion barrier. These effects will lead to a reduced fission prob- 
ability after the fission delay is over. An experimental test of 
this idea would be the observation of reduced fission probabil- 
ity or enhanced evaporation-residue survival at high excitation 
energies which cannot be explained in terms of the statistical 
model. At present there is some controversy over the need for 
fission transients. 

A number of theoretical studies iQl reproduce experimen- 
tal fission probabilities and pre-scission neutron multiplici- 
ties with transient fission widths when the viscosity increases 
with the mass of the compound nucleus. Transient fission has 
also been invoked to explain the relatively large number of 
evaporation residues measured for the very fissile ^^^Th com- 
pound nuclei formed in •^^Sh-^'^^W reactions, as compared to 
a statistical-model prediction [|10ll . Alternatively other studies 
have reproduced fission probabihties with no transient effects 
lUin and Lestone and McCalla III2I1 suggest that fission tran- 
sients are unimportant whenever the nuclear temperature is 
less than the fission barrier. Similarly, in very high excitation- 
energy data obtained in 2.5-GeV proton-induced reactions, no 
transients were needed in reproducing the measured fission 
yields IIbII . 

This paper discusses the application of the GEMINI-H- de- 
excitation model III4I1 to the description of fission in fusion 
and spallation reactions. In the latter case, the description of 
the entrance channel is provided by a coupling to the Liege 
Intranuclear Cascade model (INCL) |[l5t]. Both INCL and 
GEMINI-H- are among the most sophisticated models in their 
own fields. The present work also represents the first thorough 
discussion of their coupling. 

We compare the predictions of the models with experimen- 
tal residue yields in spallation studies and with fission and 
evaporation-residue excitation functions measured in heavy- 
ion induced fusion reactions. The choice of the observables 
was motivated by considerations about their sensitivity to fis- 
sion and by the availability of experimental data. An ex- 
haustive discussion should of course take into account other 
observables (e.g. double-differential particle spectra) and the 
competition of fission with the other de-excitation channels, 
but this is outside the scope of the present paper Therefore, 
we describe how the parameters of the statistical-decay model 
have been consistently adjusted to reproduce the data and dis- 
cuss to what extent a successful unified description of these re- 
actions has been achieved. Finally, we explore whether all the 
data can be described within the statistical model or whether 



transient fission decay widths are needed. 



II. THE MODELS 

We shall now turn to the description of the most important 
features of the models we have considered. The codes will 
not be analyzed in detail, but only the most important features 
will be outlined. 



A. GEMINI++ 

GEMINI-H- is an improved version of the GEMINI sta- 
tistical decay model, developed by R. J. Charity lllql with the 
goal of describing complex-fragment formation in heavy-ion 
fusion experiments. The de-excitation of the compound nu- 
cleus proceeds through a sequence of binary decays until par- 
ticle emission becomes energetically forbidden or improbable 
due to competition with gamma-ray emission. 

Since compound nuclei created in fusion reactions are typ- 
ically characterized by large intrinsic angular momenta, the 
GEMINI and GEMINI-H- models explicitly consider the in- 
fluence of spin and orbital angular momentum on particle 
emission. Moreover, GEMINI/GEMINI-H- do not restrict 
binary-decay modes to nucleon and light-nucleus evaporation, 
which are the dominant decay channels, but allow the decay- 
ing nucleus to emit a fragment of any mass. The introduction 
of a generic binary-decay mode is necessary for the descrip- 
tion of complex-fragment formation and is one of the features 
that set GEMINI/GEMINI-H- apart from most of the other 
de-excitation models. 

Emission of nucleons and light nuclei {Z < 2, 3 or 4, 
depending on the user's choice) is described by the Hauser- 
Feshbach evaporation formalism iflTll . which explicitly treats 
and conserves angular momentum. The production of heavier 
fragments is described by Moretto's binary-decay formalism 
yD . However for symmetric divisions of heavy compound nu- 
clei, the Moretto formalism employing Sierk's Finite-Range 
calculations lUSl [l9ll fails to reproduce the mass distribution 
of decay products (Sec. IIIICI l. However for light systems, 
the Moretto formalism works quite well Ilia, USl I20I1 and is 
still used in GEMINI-H-. Also for the heavier systems, but 
for mass-asymmetries outside of the symmetric fission peak, 
the Moretto formalism is still used. Otherwise, the total fis- 
sion yield is obtained from the Bohr- Wheeler formalism llll 
and the width of the fission-fragment mass distribution is 
taken from systematics compiled by Rusanov et al. 12111 (see 
Sec. lmCl l. 

Table U summarizes the de-excitation mechanisms featured 
by GEMINI-H-. 

The parameters of the model associated with evaporation 
have been adjusted to reproduce data from heavy-ion-induced 
fusion reactions. This is described in Ref. ISIJ in more detail, 
but we briefly list the important adjustments for that work . In 
order to fit experimental light-particle kinetic-energy spectra, 
the transmission coefficients in the Hauser-Feshbach formal- 
ism were calculated for a distribution of Coulomb barriers as- 



Table I. List of de-excitation processes featured by the GEMINI++ 
model. The symbol Z represents the charge number of the emitted 
particle and Zswitch can be chosen to be 2, 3 or 4. 



Process 


Model 


Notes 


light-particle 
evaporation 


Hauser-Feshbach fl^l 


Z ^ ^switch 


binary decay 


Moretto d 


Z > ^switch 


fission 


Bohr- Wheeler Q|] 


only in heavy 
systems 


partition in fission 


Rusanov et al. pil) 





sociated with thermal fluctuations. The nature of fluctuations 
is not entirely clear, they may be fluctuations in compound- 
nucleus shape and/or its density and/or its surface diffuseness. 
Level densities were calculated with the Fermi-gas form: 



p{E*, J) - exp hyOuyu 



(1) 



where E* is the total excitation energy, J is the spin and U 
is the thermal excitation energy after the pairing, rotational, 
and deformation energies are removed. The level-density pa- 
rameter used should be considered an effective value as no 
collective-enhancement factors are used in the level density 
formula of Eq. ([T). 

The level-density parameter a{U) is excitation-energy de- 
pendent with an initial fast dependence due to the washing out 
of shell effects following Ref. 122 and a slower dependence 
needed to fit the evaporation spectra. The shell-smoothed 
level-density parameter was assumed to have the form 



a{U) 



A 



koo - (fcoo - fco)exp (- 



-feo A 



(2) 



which varies from A/ko at low excitation energies to A/k^o 
at large values. Here fco=7.3 MeV, consistent with neutron- 
resonance counting data at excitation-energies near the neu- 
tron separation energy, and koo=l2 MeV. The parameter k de- 
fines the rate of change of a with energy and it is essentially 
zero for light nuclei (i.e. a constant a value) and increases 
roughly exponentially with A for heavier nuclei. Although 
we expect a decrease in the level-density parameter with U 
due to decreasing importance of long-range correlations with 
increasing excitation energy (due to washing out of collective- 
enhancement factors and also the reduction of the intrinsic 
level-density parameter), the strong mass dependence cannot 
be explained at present. 

The strong excitation-energy dependence of a for heavy nu- 
clei leads to increased nuclear temperatures which enhance 
very weak decay channels. For very fissile systems, these 
weak decay channels include n, p, and a evaporation and 
thus GEMINI-H- calculations predict enhanced evaporation- 
residue production consistent with some experimental data. 
These enhanced evaporation-residue yields had previously 
been interpreted as a consequence of transient fission lllOll . 
Clearly the excitation-energy dependence of a is very impor- 
tant in understanding the role of transient fission. However 



we note that for fission, the dominant decay mode in fissile 
nuclei, the yield is decreased relatively little by the increased 
temperature. 



B. INCL4.5 

The INCL model ITsIl can be applied to collisions between 
nuclei and pions, nucleons or light nuclei of energy lower than 
a few GeV. The particle-nucleus collision is modelled as a se- 
quence of binary collisions among the particles present in the 
system; particles that are unstable over the time scale of the 
collision, notably A resonances, are allowed to decay. The 
nucleus is represented by a square potential well whose ra- 
dius depends on the nucleon momentum; thus, nucleons move 
on straight lines until they undergo a collision with another 
nucleon or until they reach the surface, where they escape if 
their total energy is positive and they manage to penetrate the 
Coulomb barrier. 

The latest version of the INCL model (INCL4.5) includes, 
among other things, isospin- and energy-dependent nucleon 
potentials, an isospin-dependent pion potential and a new 
dynamical coalescence algorithm for the production of light 
clusters (up to A = 8 with the default program options). A 
comprehensive description will be published in the near fu- 
ture iH. 

The INCL model simulates a complete cascade event, its 
output being the velocities of all the emitted particles. The 
characteristics of the remnant (its mass, charge, momentum, 
excitation energy and intrinsic angular momentum) are de- 
rived from the application of conservation laws and are passed 
to the chosen de-excitation code; the latter simulates the de- 
cay of the remnant into a nuclear-stable residue plus a number 
of nucleons, nuclei and/or gamma rays. 

The INCL4.5 model is not to be considered as an ad- 
justable model. It does contain parameters, but they are either 
taken from known phenomenology (such as the matter density 
radius of the nuclei) or have been adjusted once for all (such as 
the parameters of the Pauli blocking or those who determine 
the coalescence module for the production of the light charged 
clusters). Therefore adjusting INCL4.5/GEMINI-H- on the 
experimental data basically amounts to the adjustment of 
the GEMINIh-h- parameters. One should keep in mind that 
INCL4.5 brings in its own physics features and limitations. 
For our purpose here, they essentially determine the distribu- 
tions of the remnant properties. These quantities cannot be 
compared directly with experimental data, but the predictions 
of INCL4.5 concerning those observables that can be con- 
fronted directly to experiment, namely the high energy parts 
of particle spectra, are of rather good quaUty, as it was shown 
recently 1 2411 . 



III. ADJUSTMENT OF FISSION YIELDS 

The assumption of thermal equilibrium implied by the 
statistical-decay hypothesis implies that the excited nucleus 
cannot keep any memory of the entrance channel. One of 
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Figure 1. (Coloronline)Comparisonof the distributions of excitation 
energy and spin populated in tlie ^^F+^*^Ta— ^■^""Pb fusion reaction 
for E* = 90, 150 MeV (horizontal lines) with the INCL4.5 predic- 
tion for the l-GeVp+^''*Pb spallation reaction (contours, logarithmi- 
cally spaced). The inset shows the distribution of fissioning remnants 
for the spallation reaction (same contour levels). The dashed line is 
the macroscopic yrast line from Sierk I25ll . 



the main aspirations of the GEMINIh-h- development is to 
provide a unified and coherent description of nuclear de- 
excitation in spallation and fusion reactions at the same time. 

The degrees of freedom in the model induce different char- 
acteristic dependencies of the fission width on the remnant 
spin and excitation energy, because fission is at the very least 
sensitive to spin, the fade-out of shell and collective effects, 
level densities and fission barriers. Since variations in some 
of the free parameters can produce similar effects, it is diffi- 
cult to disentangle the various contributions and to put strin- 
gent constraints on the de-excitation model just by looking 
at experiments of a single type. However, fusion and spal- 
lation reactions populate different regions of the compound- 
nucleus spin/excitation-energy plane. A comparison of the 
populations in the E*-J plane is shown in Fig. [T]for the 1- 
GeV /7H-208pb spallation reaction and the i^FH-i^iTa-^^oop^ 
fusion reaction. The spallation population is represented by 
the contours and two examples of the fusion distributions for 
E*=9Q, 150 MeV are shown by the thick horizontal shapes, 
the thickness of which is proportional to the population. The 
inset shows the distributions of remnants leading to fission in 
the case of spallation; the contour levels are the same as in the 
main plot. As a guide for ■^'^"Pb, the macroscopic yrast line 
from Sierk |25] is indicated. 

For the spallation reaction, the INCL4.5 model predicts 
average values of about 167 MeV and 16.5 h, but both dis- 
tributions are quite broad and extend up to ^ 650 MeV and 
~ 50 h, respectively. On the other hand, the fusion reactions 
we considered are characterized by higher spins and lower ex- 



citation energies. We concentrate only on complete-fusion re- 
actions, where the excitation energy of the compound nucleus 
is defined entirely from energy conservation. The require- 
ment of complete fusion restricts us to projectile bombard- 
ing energies of less then 10 MeV/A, where incomplete fusion 
and pre-equilibrium process are small. We can explore some- 
what higher excitation energies with more symmetric entrance 
channels, but high spins will still be populated. Thus it is clear 
that the comparison between spallation and fusion data repre- 
sents a promising tool to extend the predictive power of the 
model over a wide region of mass, energy and spin. 

The outlined simultaneous-fitting strategy would be mean- 
ingless if the two types of data probed separate model- 
parameter subsets; separate fits would then be just as effec- 
tive as a combined fit. This is not the case, as Fig. [T] sug- 
gests. Spallation is characterised by a large spread in the 
compound-nucleus excitation energy; however, the excitation- 
energy dependence of fission can be tested in fusion reactions, 
too, using excitation curves. Thus, both types of data sets are 
probably sensitive to the model parameters connected with the 
excitation-energy dependence of fission. A similar argument 
can be produced for the model parameters connected with the 
spin dependence. The model parameters that are sensitive to 
the mass and charge of the fissioning nucleus are probed by 
spallation through the distribution of cascade remnants, and 
by fusion when one considers neighbouring compound nu- 
clei. We then conclude that the fusion and spallation probe 
different, but overlapping areas of the model-parameter space. 
Under these conditions, simultaneous fitting is most likely ad- 
vantageous. 

The fusion reactions used in this study are listed in Table HIl 
with the range of excitation energies probed and the appro- 
priate references for the data. In most cases we have selected 
data where both the evaporation-residue and fission cross sec- 
tions have been determined. The sum of these two quantities 
gives the total fusion cross section and this is used to constrain 
the compound-nucleus spin distribution. We assume the spin 
distribution has a roughly triangular shape characterized by a 
maximum value Jq with a smooth cutoff characterized by the 
parameter SJ, i.e. 



afus{J)^TTX\2J+l) 
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The parameter Jq is determined from the total fusion cross 
section 



Cfus = 2^ O'fus(J) 



J=0 

and SJ is set to values from 3 to 10 /i, with the larger val- 
ues associated with the heavier projectiles. These values are 
similar to estimates obtained in Refs. 026142811 

Generally the fission cross section is only sensitive to the 
value of SJ at excitation energies where the fission probabil- 
ity is small and rises rapidly with J ll26ll . The values of SJ 
assumed in the following calculations are also listed in Ta- 
ble |II] 



Table II. Experimental fission and evaporation-residue data used in 
this work. 
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reaction 


E* range 
[MeV] 
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CTfus 


<5J 
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156Er 


e^Ni + «OZr 


25-82 


im 


[30] 


10 


i58Dy 


is-F+i^SLa 


68-94 


[26] 


[26] 


2.3 


leoyb 


60Ni + i00Mo 


95-249 


[SIJ 




10 


lesyb 


i«0 + ^^"Sm 


63-126 


[31,^2] 


[31,^2] 


4 


178^ 


19F+ i^^W 


54-95 


[31J 


[31J 


4.2 


188pj 


i9F + i69Tm 


50-91 


[SIJ 


[SIJ 


4.4 


193^1 


2«Si+^«^Ho 


65-249 


[Ml 




10 


200pb 


19p+ 181^3 


49-153 


[34. ia 


[34. ia 


4.7 


200pb 


^°Si+^™Er 


48-84 


[34J 


[34J 


10 


216tJj 


32s ^184^ 


125-203 


[10] 


[361 


10 


2i«Ra 


i^F + i^^Au 


39-86 


\3T\ 


[37] 


3 


2i«Ra 


30si + 186^ 


39-83 


[37] 


[37] 


10 


224^^ 


16o + 208pj, 


26-83 


[3^ 


[39,40] 
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A. Modifications of the Fission Width 



The Bohr- Wheeler fission width. 
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is sensitive to the choice of the fission barrier B and to the 
level-density parameters at and a,, associated with the saddle- 
point and ground-state configurations. The U dependence of 
the level-density parameter was initially assumed to be identi- 
cal for the ground state and the saddle point, and it is described 
by Eq. (|2]). However, the saddle-point level-density parameter 
ttf was scaled by a constant factor with respect to the corre- 
sponding ground-state level-density parameter an, to account 
for the increased surface area of the saddle-point configura- 
tion IHTr . In what follows, we will refer to the scaling factor 
as "the af/fln ratio", for simplicity. 

A number of modifications to the Bohr- Wheeler width have 
been proposed. In a one-dimension derivation of the escape 
rate over a parabolic barrier for high viscosity, Kramers ll48ll 
obtained 



For the ^^Si-n^^^Ho and ^°NiH-^™]VIo reactions, only evap- 
oration residue data has been measured. However for these 
reaction, the Jo values associated with fusion-like reactions 
are very large and the higher J values all go into fission. The 
evaporation residue yield is therefore not sensitive to Jo and 
is entirely determined by fission competition at the lower J 
values. Blann et al. f4]| termed this a saturation analysis as 
the higher J values are saturated by fission. 

For spallation reactions, we focused our efforts on proton- 
induced fission reactions on ^^^Au fill, ^''^Pb HI and ^^^V 
11441] at 1 GeV, measured in inverse kinematics with the FRag- 
ment Separator (FRS) at SIS, GSI, Darmstadt, Germany. An 
additional experimental data-set for p+^'^^Ph at 500 MeV ex- 
ists 11451] . but new measurements seem to indicate that the fis- 
sion cross section was overestimated by about a factor of 2 
mm . We decided to normalize Fernandez et al.'s total fission 
cross section to the cross section measured by the CHARMS 
collaboration, assuming that the fission distribution had been 
correctly measured. The reader should nevertheless keep in 
mind the normalization uncertainties associated with this data 
set. 

The simultaneous-fitting strategy, however, cannot be ap- 
plied to pre-scission neutron multiplicities, since there are no 
such data (to our knowledge) for spallation reactions. More- 
over, accurate modelling of the pre- and post-scission neutron 
data is likely to introduce new ingredients and parameters re- 
lated to the the saddle-to-scission descent (e.g. the viscosity 
of the saddle-to-scission motion). It is not obvious whether 
the extra constraint provided by the new data would overcom- 
pensate the increase in the number of degree of freedom and 
lead to an effective decrease of the model uncertainty, all the 
most so in the lack of relevant spallation data. Thus, we will 
not discuss pre-scission neutron multiplicities in this work. 
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where 7 is the magnitude of the viscosity, uj is frequency asso- 
ciated with the inverted parabolic barrier, and the factor scal- 
ing the Bohr- Wheeler decay width is less than unity. Now uj is 
not expected to be a strong function of mass or spin, and if 7 is 
also constant, then the Kramers and the Bohr- Wheeler values 
differ by approximately a constant scaling factor. For this rea- 
son we have allowed a constant scaling to the Bohr- Wheeler 
width. 

Lestone 11491] developed a treatment of fission which explic- 
itly included the tilting collective degree of freedom at saddle 
point. Tilting is where the compound nucleus's spin is not 
perpendicular to the symmetry axis. For strongly-deformed 
objects like the saddle-point configuration, this costs energy 
and thus decreases the fission probability. The decay width 
becomes 
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where the summation is over K, the projection of the spin on 
the symmetry axis and 
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and /|| and I± are the saddle-point moments of inertia parallel 
and perpendicular to the symmetry axis, respectively. In this 
work, the moments of inertia as well as the spin-dependent 
saddle-point energies were taken from the finite-range calcu- 
lations of Sierk 11251] . Deviations from the Bohr- Wheeler value 
are largest for the highest spins and thus the Lestone modifi- 
cation will be more important in fusion reactions. 

We have tried to reproduce simultaneously fission cross 
sections from fusion and spallation experiments by: 
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Figure 2. (Color online) Comparison of GEMINI++ predictions to 
the experimental evaporation-residue and fission excitation functions 
for the ^''F+i^^Ta reaction. 



• Adding a constant to the Sierk fission barriers for all 
spins. 

• Scaling the decay width by a constant factor. 

• Adjusting the af/on ratio. 

• Using either the Bohr- Wheeler or the Lestone formal- 



• Introducing a constant fission delay. 



B. Fission Probability 

Examples of fits to the i^F-n^^^Ta^^^Pb fission and 
evaporation-residue excitation functions are shown in Fig.|2l 
As the sum of these quantities (the fusion cross section) 
is fixed in the calculations, the degree to which the fission 
probability is reproduced is best gauged by the fit to the 
smaller quantity, i.e. o-fis at low bombarding energies and 
(Ter at the higher values. Good fits were obtained with 
Fbw X 2.46, Of/ On ~ 1.00 (long-dashed curves), Fbw x 
1.00,af/an = 1.036 (solid curves), FLestone x 7.38,af/an = 
1.00 (dotted curves), and FLestone x 1.00,af/an = 1.057 
(short-dashed curves). The Fbw x 2.46,af/an = 1.00 calcu- 
lation is also almost identical to a Fbw x 1.00, af/cn = 1.00 
calculation (not plotted) obtained with the Sierk fission barrier 
reduced by 1 .0 MeV. With an even larger barrier reduction fac- 
tor, one could arrive at a solution where the decay-width scal- 
ing factor is less than unity and consistent with the Kramers' 
scaling factor in Eq. (O. 

As it is impossible to distinguish these different ways of 
modifying the fission probability from the fusion data alone, 
we now consider the constraint of adding the spallation data 
to the analysis. In Fig. |3] we show the equivalent calcula- 
tions for the mass distributions of the products of the 1-GeV 
pH-^^^Pb spallation reaction. Of all these possibilities, the 
Fbw X 1.00,af/an = 1.036 calculation reproduces the yield 
of the fission peak best. This highlights the significant re- 
duction in fitting-parameter ambiguity that can be obtained 
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Figure 3. (Color online) Comparison of experimental and calculated 
residual mass distributions for the 1-GeV p-F^"*Pb reaction. Pre- 
dictions of the INCL4.5-GEMINI++ code are shown for different 
adjustments of the fission width. Experimental data from Ref. 14311 . 



by simultaneously fitting heavy-ion and spallation data. We 
also note that the Fbw x 2.46, Of/on = 1.00 calculation 
and the reduced fission barrier calculation with Fbw x 1.00, 
af/fln = 1-00 (not shown) were again identical. Thus, while 
many of the ambiguities in the fitting-parameters have been 
removed, the ambiguity between the effect of the magnitude 
of the fission barrier and the Kramers scaling factor was not 
not lifted with simultaneous fitting. We will continue using the 
Fbw X 1.00, a^ja^ ~ 1.036 calculation as our best fit to both 
sets of experimental data, but it should be noted that with re- 
duced fission barriers, an equivalent solution with a Kramers 
scaling factor (< 1) can also be obtained. We are just not 
able to constrain the magnitude of the Kramers factor from all 
these data. 

Comparison of GEMINIh-h- predictions to the experimen- 
tal fission and evaporation-residue excitations functions listed 
for the heavy-ion-induced fusion reactions in Table |II] are 
shown in Figs. |4] to |2l The solid curves show the predic- 
tions with af/an=1.036, Sierk fission barriers, and no scal- 
ing of the Bohr- Wheeler decay width. For spallation. Figs. [8] 
to [To] show the comparison between measured and calculated 
residue mass distributions. Finally, Fig.[TT]shows the compar- 
ison between measured IdOII and calculated excitation curves 
for the fission cross section in proton collisions with ^^^Ta, a 
low-fissility target. 

The central result is that it is possible to reproduce the to- 
tal fission cross section for all the studied spallation reactions 
by adjusting only one free parameter, namely the af/fln ratio, 
which was set equal to 1.036 in our calculations, whilst the 
global scaling of the fission width and of the fission barrier 
were kept equal to 1; no Lestone coiTection was introduced. 
A global scaling of the fission width is roughly equivalent 
to a reduction of the barrier height, but in both cases, these 
adjustment alone do not fit the data. The adjustment of the 
af/fln ratio, on the other hand, is characterized by a different 
excitation-energy dependence, which is better suited for the 
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Figure 4. (Color online) Comparison of experimental and calculated 
fission and evaporation-residue excitation functions for the indicated 
reactions. Solid lines: Bohr- Wheeler fission width, at/an — 1.036, 
no fission delay. Dashed lines: Lestone fission width, at/an = 
1.057, 1-zs fission delay. Dotted lines: Bohr-Wheeler fission width, 
energy-dependent effective af/a„ ratio with r= 1.0747. 



description of fission from spallation remnants. The Lestone 
correction, which suppresses the fission width at high spin, 
does not have a large effect on spallation data, some 80% at 
most. This is due to the small angular momenta generated in 
the intranuclear cascade. 

For the heavy-ion-induced fusion data, the GEMINIh-h- 
predictions are also generally quite good, however, there are 
a couple of reactions where significant deviations are found. 
Firstly, for the i^On-isoSm ->i68Yb and iSFH-i^SLa-^iss^y 
reactions shown in Figs.|6jb) and|6jc), the fission cross sec- 
tion is overestimate by almost an order of magnitude. In com- 
parison, the Ni-induced reactions in Figs.|7jb) and|7Jc) mak- 
ing similar mass compound nuclei are reproduced much bet- 
ter. It is difficult to understand these ^^O- and ^^F-induced 
reactions and previous attempts also failed to reproduce the 
data 12611. For instance at the highest bombarding energies, 
the excitation energies probed in these reactions overlap those 
in the Ni-induced reactions. Similarly the predicted angular- 
momentum region over which the fission yield is determined 
in the O- and F-induced reactions is similar to that in which 
the residue yield is determined in the Ni-induced reactions. 
Unless there are significant non-fusion processes, such as pre- 
equilibrium or incomplete fusion occurring, these data suggest 
an entrance-channel dependence of the fission decay proba- 



Figure 5. (Color online) Same as Fig.|4] 
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Figure 6. (Color online) Same as Fig.|4] 
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Figure 7. (Color online) Same as Fig.|4l 
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Figure 8. (Color online) Residue-mass distribution for the p+ Au 
reaction at 1 GeV. Predictions of the INCL4.5-GEMINI++code are 
shown for different adjustments of the fission width. Experimental 
data from Ref. HI. 



bility which would violate the compound-nucleus hypothesis. 
However we see no evidence of such an effect for the O- and 
F-induced induced reactions with heavier targets. Clearly our 
understanding of fission for A < 170 is lacking and more 
studies are needed. 

The other case where the GEMINI-n- predictions fail is 
for the ''^SH-iS'^W-^^ieT^'jj reaction in Fig. [Sja). Here the 



< 



10 



rBW'af/an=''-036 



rLestone.af/an="'-057,t,,„y=1ZS 



rB^y,r=1.0747 

'^Lestone' *clelay~^ ^®' •"- ' •' 





40 60 80 100 120 140 160 180 200 



Figure 9. Same as Fig.[8]forp+^°*Pb at 500 MeV. Experimental data 
from Refs. 14a. Dill . The experimental fission cross sections have 
been multiplied by a factor of 146/232 (see text for details). 
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Figure 10. (Color online) Same as Fig.[l]forp+^^*U at 1 GeV. Ex- 
perimental data from Refs. I44l . l52l . 15311 . 



evaporation-residue cross section is exceedingly small (~ 
0. 1 mb) and is overpredicted by almost an order of magnitude. 
However, the calculations gets the excitation-energy depen- 
dence of the cross sections correct which previous calculations 
could not do without invoking an excitation-energy depen- 
dence of the dissipation strength jlOtl . In our calculations the 
predicted excitation-energy dependence of the residue cross 
section is a consequence of the assumed excitation-energy de- 
pendence of the level-density parameter |5[|. Low probability 
events in the statistical model are generally quite sensitive to 
the statistical-model parameters. In this case, it was demon- 
strated that the residue yield is very sensitive the absolute 
value of the level-density parameter and its excitation-energy 
dependence l^. For example, the residue yield is increased 
by 2-3 orders of magnitude when the level-density parameter 
is changed from a = A/7.3 to A/11 MeV^^. Further re- 
finement of the value of this parameter at the larger excitation 
energies probed in this more symmetric fusion reaction may 
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Figure 11. (Color online) Excitation curve for the fission cross sec- 
tion inp+^**^Ta. Predictions of the INCL4.5-GEM1N1++ code are 
shown for different adjustments of the fission width. Experimental 
data from Ref. Isoll. 



be needed in the future. 

Alternatively, there is evidence that for this mass region, 
quasi-fission completes with fusion reactions even at the 
lower i waves associated with evaporation residue produc- 
tion OTI 15411 . In the case of ^^^Ra compound nucleus in 
Figs. |5lb) and|5lc), Berriman et al. have indicated that both 
the ^^Fh-^^^Au and "^"Si-n^^'^W reactions have reduced evap- 
oration residue cross sections due to quasi-fission competi- 
tion IBTT . Therefore this higher mass region for the heavy- 
ion reactions is subject to more uncertainty in constraining 
the statistical-model parameters. 



C. Fission-Fragment Mass Distributions 

Previous treatments of the fission-fragment mass distribu- 
tion have assumed thermal models where the mass division 
is determined at either the saddle-point (Moretto's formal- 
ism) or at the scission configuration 15511 . In reality, a com- 
plete description probably requires one to follow the trajectory 
from saddle to scission including fluctuations, for example by 
Langevin simulations iDapTll . However for large-scale sim- 
ulations, this is too time consuming, so a simpler and faster 
procedure for determining mass division is required. 

Experimental mass distributions for heavy-ion-induced 
fusion-fission reactions are shown in Figs. [12] and [13] In these 
figures, the experimental fission-fragment masses was not di- 
rectly measured, but rather the ratio of primary masses (be- 
fore post-scission particle evaporation) was inferred from ei- 
ther the ratio of the measured fission-fragment velocities or 
kinetic energies. The absolute primary mass was assumed to 
be equal to the compound-nucleus mass, which of course ig- 
nores the pre-scission evaporation of light particles. However, 
the distributions simulated by GEMINIh-h- (the curves in the 
figures) were analyzed in the same manner as the experimen- 
tal data and thus contain the same deficiencies. 




primary 



Figure 12. (Color online) Comparison of experimental (data points) 
and simulated (curves) distributions of the primary fission-fragment 
masses for the ^®*Pb compound nucleus formed in the ^"^O+^^^W 
reaction at the two indicated excitation energies. The total fission 
cross section has been calculated using the Bohr- Wheeler formal- 
ism, ftf/fln = 1.036, without any fission delay. The dotted curves 
was obtained using the Moretto formalism with Sierk's conditional 
barriers to define the mass distributions. The solid and dashed curves 
were obtained using the Rusanov systematics using the saddle and 
scission-point temperatures, respectively. 
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Figure 13. (Color online) Comparison of experimental (data points) 
and simulated (curves) distributions of the primary fission-fragment 
masses for the ^^''Ra compound nucleus at E* ~ 60 MeV formed 
in the ^^C+^""'Pb, ^^F+^'-*^Au and ^"Si+^^^^W reactions. The total 
fission cross section has been calculated using the Bohr- Wheeler for- 
malism, af/a„ = 1.036, without any fission delay. To aid in viewing, 
the data and curves have been scaled by the indicated amounts. 



The data in Fig.[l2]is for the iSOH-i*2w->i98pb reactions. 
The relative mass distribution were obtained from Ref. 15811 
and absolute normalization was achieved using the fission 
cross sections measured in Ref. I159ll . The data in Fig. [13] are 
for the •^^^Ra compound nucleus at E* r^ 60 MeV with three 
entrance channels; i2cH-204pb^ isFh-I^^Au, and ^Ssin-isew. 
The fission excitation function for the latter two are shown in 
Figs.|5]:b)and|5];c). 

In Fig. [12] the dotted curves show the mass distribution de- 
termined from Moretto's formalism using interpolated values 
of Sierk's finite-range calculations for the conditional barriers 
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Ill4tl . The total fission width (the total width for all mass di- 
visions associated with the peak in the mass distribution) was 
normalized to the Bohr- Wheeler value for these calculations. 
Therefore in this figure, only the shape of the mass distribu- 
tion is determined from the Moretto formalism. Clearly these 
distributions are much wider than the experimental quantities. 
This is quite typical of other cases where a peak exist in the 
mass distribution at symmetry. See for example the study of 
^^^Eu compound nuclei in Ref. 1I6OI1 and the light, but high 
spin ^^"Sn compound nuclei studied in Ref. 1I6II1 . For less 
fissile nuclei where the mass distribution has a minimum at 
symmetric division, the Moretto formalism (with Sierk barri- 
ers) give a much better description of the experimental data. 
See for example the studies of i"ln lO, ^°^Rh and ^o^Ag 
HI] and "^^Br Q compound nuclei. 

The cause of this inadequacy for the heavier systems could 
either be an incorrect asymmetry dependence of Sierk's con- 
ditional barriers or a failure of Moretto's formalism. The latter 
predicts the asymmetry distribution at the ridge line of condi- 
tional saddle points and assumes that the mass asymmetry is 
unchanged during the descent from saddle to scission. 

As an alternative to using Moretto's formalism, we have 
used the systematics of fission-fragment mass distributions 
complied by Rusanov et al. 112111 . The mass distribution is as- 
sumed to be Gaussian and its variance is parameterized as 
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where 



A^V 



j^2 is the second derivative of the potential en- 
ergy surface with respect to the mass-asymmetry deformation 
parameter(7^ ~ 2 ^~^ where Ai and A2 is the mass divi- 
sion). This quantity is parameterized as a function of the fis- 
sility Z'^ /A and spin J. The quantity T is the nuclear temper- 
ature where 
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T 



dlnp 
dU ■ 



Rusanov et al. considered three parameterizations of ^-^, 
with three different temperatures. Either 1) the temperature 
of the fission nucleus at the saddle-point is used, but no pre- 
saddle light-particle evaporations are allowed, 2) as above, but 
pre-saddle evaporations are allowed, or 3) the temperature at 
the scission point is used. The first of these is not realistic 
and was not considered and the second is basically consistent 
with the ideas of the Moretto formalism. The latter two can 
be called saddle-point and scission-point models where the 
mass distributions are both determined thermally. In these two 
cases, the quantity ^^^ should be identified with the asymme- 
try dependence of the potential-energy surface at the saddle 
and scission points, respectively. 

In GEM1N1-H-, once fission is decided for an event, evapo- 
ration during the saddle-to-scission transition is allowed. This 
is important for the scission model, as we need to determine 
the temperature at the scission point. The saddle-to-scission 
evaporation is treated in a simplified manner using spherical 
level densities and transmission coefficients in the Weisskopf- 
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Figure 14. (Color online) Comparison of experimental (data points) 
and simulated (curves) distributions of the fission-fragment masses 
for the 1-GeV p+Ph spallation reaction. The blue dotted curve was 
obtained using the Moretto formalism with Sierk's conditional barri- 
ers to define the mass distributions. The black solid and red dashed 
curves were obtained using the Rusanov systematics using the saddle 
and scission-point temperatures, respectively. 



Ewing evaporation formalism with the deformation-plus- 
rotational energy removed from the total excitation energy. 
The deformation-plus-rotational energy of the scission con- 
figuration is determined as the sum of fission-fragment kinetic 
energy from Viola's systematics ll62ll and the fission Q-value. 
Evaporation during the saddle-to-scission transition occurs for 
a period proportional to the difference in energy between the 
saddle and scission points, i.e. 



U. 



"'Ssl^-^St 



E,, 



d 



consistent with large viscosity. The parameter fcss is related to 
the magnitude of this viscosity was fixed to fcss=l zs/MeV by 
fitting pre-scission neutron multiplicities from Ref. fT]. 

The solid and dashed curves in Fig.[T2]show the predictions 
with the Rusanov saddle and scission-point systematics, re- 
spectively. These predictions are almost identical, and for the 
lowest excitation energy, the curves completely overlap and 
cannot be distinguished. This is not surprising as both Ru- 
sanov systematics are fits to "^He induced fission and fusion- 
fission data including the data set of Fig. [12] However, at 
higher excitation energies such as those sampled in spalla- 
tion reactions, the two systematics give quite different pre- 
dictions as the thermal excitation at scission increases much 
more slowly with compound-nucleus excitation than does the 
saddle-point value [Tj]. Figure [T4l compares, the two system- 
atics for the 1-GeV pH-Pb spallation reaction. In this case, the 
predicted mass distribution obtained with the scission system- 
atics (dashed curve) is too narrow, while the saddle systemat- 
ics (solid curve) gives good agreement. 

The success of Rusanov's saddle-point systematics thus 
suggests that the fission mass division is determined quite 
close to the saddle-point configuration. It addition it indicates 
that the Moretto formalism is still applicable for near symmet- 
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Figure 15. (Color online) Comparison of asymmetry dependences of 
conditional barriers for ^''^Tb and ^^*Hg nucleus at J=0. The solid 
curves are the predictions from Sierk's finite-range calculations. The 
dashed curves are parabolic functions with curvatures taken from the 
Rusanov systematics and with the symmetric fission barriers taken 
from the finite-range calculations. The mass-asymmetry coordinate 
is defined in terms of Ai and A2, the two masses following binary 
division. 



ric divisions of heavy nuclei. However it should not be used 
with Sierk's conditional barriers in this region. 

The differences between the mass-asymmetry dependence 
of Sierk's conditional barriers and the Rusanov systematics 
are shown directly in Fig.[T5]for ^^^Tb and ^^"'Hg compound 
nuclei at J=0. The dashed curves are parabolic functions with 
curvatures from the Rusanov systematics and with the sym- 
metric fission barriers from Sierk's calculations. The Rusanov 
results have larger curvatures at symmetry than Sierk's pre- 
dictions and thus give narrower fission-fragment mass distri- 
butions. The differences between Sierk's predictions and the 
Rusanov systematics is much larger for the heavier ^^*Hg nu- 
cleus. For even heavier nuclei, the asymmetry coordinate in 
the finite-range calculations becomes undefined as the saddle- 
point configuration has no well defined neck |63i|. If this is the 
case, then the Moretto formalism is no longer applicable for 
these systems and the mass asymmetry is determined during 
the descent from saddle to scission. In such cases the inter- 
pretation of the Rusanov systematics in terms of a Moretto- 
type may be suspect. We note that the Z'^/A dependence of 
d^V/dr]^ in Rusanov systematics has an abrupt slope change 
at Z'^/A=24 possibly related to this effect. However, even in 
the p+V spallation reaction we produce Z'^ /A ratios that are 
below this value. 

The Rusanov saddle systematics was used for the other 
spallation predictions in Figs. |8] to [TOl and gives quite good 
agreement. However, for the p+^'^^V reaction in Fig.[TOl the 
simulation fails to reproduce the small shoulder in the fission 
mass distribution for higher mass. The Rusanov systematics 
only gives the width of the distribution and will not predict 
finer structures linked to shell effects, such as this. 

In Fig. [13] the simulated mass distributions (from the sad- 
dle systematics) for the ^^®Ra compound nuclei reproduce the 
data reasonably well with the exception of the ^^CH-^°*Pb data 
where experimental distribution is somewhat narrower. Ber- 




Figure 16. (Color online) Same as Fig. [3] for the Bohr- Wheeler fis- 
sion width, Of/an = 1.036 and three different values of the fission 
delay. 
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Figure 17. (Color online) Same as Fig. [3] for the Lestone fission 
width, Of/fln = 1.057 and three different values of the fission delay. 



riman et al. ||37t| suggest that ^^CH-'^°''Pb data is all fusion- 
fission while the ^^Fh-^^^Au and ^°Si-i-^^^W data both con- 
tain quasi-fission contributions making the mass distributions 
wider This would imply that for the more massive compound 
nucleus, the Rusanov systematics overestimate the width of 
the statistical fission mass distributions as many of the heavy- 
ion data used in these systematics have contributions from 
quasi-fission. However, the Rusanov systematics also con- 
tains the lower-spin •^He-induced fission data in this mass re- 
gion and here quasi-fission is expected to be absent. Thus 
spallation mass distributions which sample lower spins are not 
expected to suffer from this problem. 



D. Fission Delays 

Apart from the lighter compound nuclei, we have demon- 
strated that a unified description of fission widths in fusion 
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and spallation reactions can be obtained. The solution is 
unique, apart from an ambiguity between the height of the 
fission barrier and the Kramers scaling factor However we 
will now show that another ambiguity arises when fission tran- 
sients are considered. To show the sensitivity of predictions to 
transients, we have incorporated a simple implementation of 
these in GEMINI++; the fission width is set to zero for a 
time tdeiay, after which it assumes its asymptotic value. Dur- 
ing this fission-delay period, the compound nucleus can decay 
by light-particle evaporation and intermediate-mass-fragment 
emission. The fission delay is expected on theoretical grounds 
to be logarithmically dependent on nuclear temperature 161], 
but this weak dependence (and any mass dependence) has 
been neglected in a first approximation. Figure [16] shows the 
dependence of the predicted mass distributions for the 1-GeV 
/9H-^"^Pb reaction with tdeiay=0, 1, and 10 zs. Even a small 
1 zs delay has a large effect on the yield in the fission peak. 
Therefore, the spallation reactions should be quite sensitive to 
the fission transients. Tishchenko et al. also expected large re- 
ductions is the fission probability in 2.5-GeV/7H-^^^Au, ^"^Bi, 
238y reactions due to fission transients; however, they were 
also able to reproduce the fission yield within the standard 
statistical-model framework II13I1. 



Jing et al. 116411 find the effect of increasing the fission de- 
lay can be largely counteracted by increasing the value of the 
o-i/da parameter Both parameters have little effect on the 
fission probability at low excitation energies. However with 
increasing excitation energy, the fission probability becomes 
ever more sensitive to both fjeiay and af/a„. Even with the 
large range of excitation energies explored in this work, we 
found it is impossible to break the ambiguity between tdeiay 
and af/fln. To illustrate this. Fig. [TS] compares the ^oop^ f^_ 
sion data to GEMINI-H- calculations with fission delay for 
both the Bohr-Wheeler [Fig. [Ha)] and Lestone [Fig. [Tib)] 
formalisms. The values of tdeiay and af/a„ listed in these fig- 
ures were obtained by reproducing the fission cross section in 
the 1-GeV pH-Pb spallation reaction. For the Bohr- Wheeler 
case in Fig. [T8l' a). one see that the calculations with tdeiay = 1 
and 10 zs are almost identical and within 30% of the experi- 
mental values. The calculation with fdeiay = zs fits the data 
somewhat better, but all calculations can be deemed accept- 
able. 



For the Lestone formalism in Fig.[T8jb), the inclusion of a 
delay with tdeiay > 1 zs improves the agreement with the data. 
As in Fig. [TSl a). the calculations with tdeiay > 1 zs are again 
almost identical. The Lestone prescription with fission delay 
also allows good agreement with the other data sets we have 
considered, see the dotted-curves in Figs. [4] to [7] (fusion) and 
Figs. [8] to [To] (spallation) which were obtained with at/an = 
1.057 and tdeiay=l zs. Calculations with the larger tdeiay values 
produce a similar level of agreement. It is thus clear that the 
magnitude of the fission transients cannot be deduced from 
the fission probability alone. 
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Figure 18. (Color online) Comparison of GEMINI++ predictions 
using (a) the Bohr- Wheeler and (b) the Lestone fission formalism to 
the experimental evaporation-residue and fission excitation functions 
for the ^''F+^*^Ta reaction. The curves are labeled by the tdeiay and 
af/a„ values obtained from fitting the fission cross section for the 
1-GeV p+Ph reaction. 



IV. FISSION AT VERY HIGH EXCITATION ENERGY 

Fission cross sections in fusion and spallation reactions 
are dominated by the most densely populated regions of the 
compound-nucleus E*-J plane (Fig. [1]). The successful re- 
production of these data thus indicates that the GEMINI-H- 
model gives an efficient description of fission from compound 
nuclei with excitation energy up to ^ 300 MeV and spin up 
to - 60 ft. 

It is possible to probe beyond this region if one considers 
other types of data. Tishchenko et al. |[65ll studied proton- 
induced spallation reactions at 2.5 GeV on gold, bismuth and 
uranium targets. They measured the fission probability in co- 
incidence with neutron, hydrogen and helium multiplicities, 
which can be used to reconstruct the excitation energy af- 
ter the intranuclear cascade. They were able to reproduce 
the measurements with an old version of the INCL-GEMINI 
model by tuning the value of af/on on a system-by-system 
basis, ranging from 1.000 for the uranium target to 1.022 for 
the gold target. These af/un values are smaller than those dis- 
cussed in the present work. 

Fig. [19] indeed shows that our candidate parameter sets 
(Bohr- Wheeler, af/a„ ~ 1.036; Lestone, af/a„ = 1.057, 
tdeiay — 1 zs) largely overestimate the fission probability at 
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Figure 19. (Color online) Fission probability as a function of the 
excitation energy of remnants in 2.5-GeV p+Au, Bi and U reac- 
tions. Solid lines; Bohr-Wheeler width, af/a„ — 1.036. Dashed 
lines: Lestone, at/an — 1.057, ideiay = 1 zs. Dotted line: Bohr- 
Wheeler width, af/on — 1.02. Dash-dotted line: Bohr-Wheeler 
width, af/a„ = 1.00. 



high excitation energy deduced by Tishchenko et al. for p+U. 
Note that the fission probability is well reproduced up to a 
few hundred MeV, which is coherent with the results of the 
previous section. The shape of the curve is indeed very sensi- 
tive to the value of a{/a„, as Fig. [T9l shows. We can interpret 
this result an indication of the fact that while a large value of 
Of/on value is appropriate at low excitation energy, the effec- 
tive flf/on value at high excitation energy should be smaller. 

An energy-dependent af/on ratio can naturally appear, 
among other things, as a consequence of the fade out of long- 
range correlations. To obtain a better reproduction of the the 
Tishchenko data, we have considered a simple refinement of 
the formula for the level-density parameter at saddle point, 
Eq. (|2]i, as follows: 



at ([/) 



A 



r {koo - ko) exp -/ 



-fcf, A 



(4) 



The r variable, which replaces the af/on ratio, is a free param- 
eter that describes the difference in the effect of long-range 
correlations for the saddle point. In the limit of zero excita- 
tion energy, Eq. (|4]i leads to 
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while for [/ — > oo, Of = On = A/ koo- The value of r thus de- 
termines the Of/fln ratio at low energy. We expect on physical 
grounds that r should be slightly larger than one, to reflect the 
increase in surface area and an enhanced collective enhance- 
ment of the saddle-point configuration. This would also lead 
to Of/on > 1 at small U. The parameter /, on the other hand, 
expresses the different fade-out rate of long-range correlations 
at the saddle point compared to the ground state. This quan- 
tity is essentially unconstrained by experimental data. How- 
ever, we observe that, from Sec. IIII BI the approximation of 



an energy-independent af/an ratio is a good one at low ex- 
citation energies, since we can successfully reproduce fission 
cross sections in fusion and spallation. We impose this condi- 
tion by requiring that 



diaf/an) 



dU 



= 0. 
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This introduces a correlation between the parameters / and r: 



I 



r (fcpo - fcp) 
rko 



(5) 



There is no a priori reason to expect that the fade-out rate at 
the saddle point (described by /) should be correlated with 
the Of/fln ratio at low energy. We make this assumption on a 
phenomenological basis. Note that Eq. ^ implies that / < 1 
for r > 1, i.e., that long-range correlations should fade out 
more slowly at the saddle point than in the ground state. One 
should also note that k has a very strong mass dependence 151] 
and therefore the modification of Of /«„ with excitation energy 
is much stronger for the p+U reaction compared to the lighter 
systems. One can indeed see from Fig.[T9]that this is the sys- 
tem that requires the biggest modification from our previous 
solution. 

We can finally determine the value of the r parameter by 
requiring, for example, that the fission cross section for 1- 
GeV p+^'^^Ph be correctly reproduced. For a Bohr- Wheeler 
width without fission delay, this condition yields r = 1.0747, 
which corresponds to af/an = 1.051 for t/ = 0. For a Le- 
stone width with a 1-zs fission delay, we get r = 1.1 and 
af/on = 1.069 at C/ = 0. Fusion-fission and spallation-fission 
are not severely affected by this modification, as shown by 
the dotted and dashed-dotted curves in Figs. l4UTn The re- 
sulting fission probability curves to the 2.5-GeV reactions are 
shown in Fig.|20] We have good quantitative agreement up to 
~ 400 MeV and we can qualitatively reproduce the decrease 
of fission probability with excitation energy for the uranium 
target. This proves that the fission probability at very high 
excitation energies is indeed sensitive to the fade-out of col- 
lective effects at saddle point. 

We have thus shown that one can obtain similar quality 
of agreement with or without a fission delay, provided that 
one increases the value of r. Therefore, while we agree with 
Tishchenko et al. 11311 that no transients are needed to explain 
their data, one cannot rule out the presence of fission tran- 
sients as well. However, the above calculations and conclu- 
sions should be taken with caution. Level densities and fission 
probabilities are dramatically sensitive to all their ingredients 
at such high excitation energies. We have indeed observed 
that the shape of the curve can also be modified by e.g. in- 
troducing surface terms (ex A'^^^) in the level-density param- 
eter formula, or by considering different functional forms for 
the fade-out of long-range correlations. Even the competition 
with evaporation cannot be neglected. Finally, very high exci- 
tation energies will eventually give rise to other phenomena, 
such as nuclear expansion and multifragmentation, which are 
not accounted for in our framework. With all these consider- 
ations in mind, we conclude that pursuing perfect agreement 
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Figure 20. (Color online) Same as Fig. [191 for a Bohr- Wheeler width 
(solid lines) or a Lestone width with a 1-zs fission delay (dashed 
lines) with an energy-dependent effective af/a„ ratio, Eq. lHJl. 



between calculations and measurements of fission probabili- 
ties for E* > 500 MeV is useless for our understanding of the 
physics of de-excitation. 



V. CONCLUSIONS 

We have described the first coupling of the Liege In- 
tranuclear Cascade model with the GEMlNl-n- compound- 
nucleus de-excitation model. The fission probability was 
calculated using Sierk's finite-range liquid-drop fission bar- 
riers ll25fl and the excitation-energy-dependent level-density 
parameters of Ref. |l5[] adjusted to reproduce experimental 
kinetic-energy spectra of light particles. The latter were very 
important to obtain the correct excitation-energy dependence 
for heavy systems. It was demonstrated (Sec. HIH i that it is 
possible to describe fission cross sections from spallation and 
heavy-ion fusion reactions for 160 < A < 230 within the 
same framework. Spallation and fusion reactions populate 
different regions of the compound-nucleus parameter space, 
and thus they probe different, but overlapping areas of the 
model-parameter space. Thus, the simultaneous fitting of the 
statistical-decay model parameters to spallation and fusion ac- 
tually allows one to lift some of the degeneracy. However, 
even with the large range of spin and excitation energy stud- 
ied, no unique parameter set could be obtained and there re- 
mained some ambiguities in the choice of parameters. In par- 
ticular, the effect of an increasing fission delay associated with 
fission transients could be offset by an increase in the param- 



eter flf/on, the ratio of level-density parameters at the saddle- 
point and ground-state configurations. In addition, modifica- 
tions to the height of Sierk's fission barrier could be offset by 
scaling of the fission decay width which could be associated 
with the Kramers scaling of the Bohr- Wheeler decay width 
due to friction. In spite of these ambiguities, we present two 
sets of statistical model parameters suitable for predictions of 
fission probabilities for spins up to 60 h and excitation ener- 
gies up to ^ 300 MeV. 

From the study of the width of the fragment-fragment mass 
distributions in both fusion and spallation reactions, we were 
able to differentiate between the systematics compiled by Ru- 
sanov et al. based on thermal distributions at either the saddle 
or scission-point (Sec. IIIICI ). Only the saddle-point system- 
atics provided good reproduction of the experimental data in 
both types of reactions, thus suggesting the fission mass divi- 
sion is determined close to the saddle point. The asymmetry 
dependence of the saddle-point conditional barriers in the Ru- 
sanov systematics is stronger than Sierk's prediction, which 
produces very wide fission-fragment mass distributions when 
incorporated in the Moretto formalism. Further indications 
could in principle be extracted from the study of other ob- 
servables, such as pre- and post-scission neutron multiplici- 
ties. However, this would happen at the expense of introduc- 
ing new parameters and ingredients for the description of the 
saddle-to-scission dynamics, which would be difficult to con- 
strain due to the lack of relevant data for spallation reactions. 

We have proven (Sec. |IV| ) that we can qualitatively de- 
scribe fission probabilities at excitation energies higher than 
300 MeV by accommodating different fade-out rates for the 
ground-state and the saddle-point configurations. However, 
we cannot exclude that other solutions are possible, given the 
uncontrollable sensitivity of the predictions of the model to a 
large number of its ingredients. Thus, we conclude that the 
theoretical uncertainties on fission probabilities at very high 
excitation energy are too large to permit drawing strong con- 
clusions about the physics of highly excited nuclei. 
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